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Abstract 

We simulate the Ising model on a set of fixed random 0^ graphs, which corresponds to a quenched 
coupling to 2D gravity rather than the annealed coupling that is usually considered. We investigate 
the critical exponents in such a quenched ensemble and compare them with measurements on dy- 
namical (j)'^ graphs, flat lattices and a single fixed (p^ graph. 

Submitted to Phys Lett B. 



^ Address: Sept. 1993 - 1994, Permanent Address: Maths Dept, Heriot-Watt University, Edin- 
burgh, Scotland 



It has recently become possible to calculate the critical exponents for various spin models on a par- 
ticular class of dynamical connectivity lattices using both the methods of conformal field theory ||l[ and 
matrix models M. Remarkably, it is even possible to solve the Ising model exactly in the presence of 
an external field on such a lattice, which has proven impossible on fixed 2D lattices [g|. The key to the 
calculations, whether in the continuum conformal field theory formalism or the discrete matrix model 
formalism, lies in the observation that putting the model on a dynamical lattice and allowing it to in- 
teract with the lattice is equivalent to coupling the model to 2D quantum gravity. This has the effect of 
"dressing" the conformal weights Aq of the critical continuum theory on a flat 2D lattice, to give new 
conformal weights A 

A-Ao = ~^A(A-l), (1) 

where 

1 



^^- — {V25~c-VT~c) (2) 

and c is the central charge of the theory in question. This, in turn, will effect the critical exponents of 
theory. These can be calculated from the conformal weights of the energy and spin operators A^ and Ao- 
respectively [^ which give a (the specific heat exponent) and (3 (the magnetization exponent) 

1-2A, 
a 



1-A, 



in the usual fashion. Given a and (3 from the conformal field theory we can then use the various standard 
scaling relations to calculate the full array of exponents. For the exactly soluble case of the Ising model 
on a dynamical lattice the standard scaling relations can be shown to hold, so it is not too great a leap 
of faith to assume their validity for other models. 

There is a growing body of numerical evidence confirming the new dynamical exponents in simulations 
l^, H, 0, both on dynamical triangulations and their dual dynamical <j)^ graphs with both spherical and 
toroidal topology lattices [^ The matrix model approach suggests that any dynamical polygonalization 
of a surface should give the same exponents for a given model, so there should be universality in this 
sense. If we think of the numerical simulations purely as an exercise in statistical mechanics it is clear 
they represent taking an annealed average over the various random graphs on which the spin models live. 
The partition function of the Ising model on dynamical lattices, for instance, is 

Z = ^^exp/3 I ^ G,,a,a, +i7^a, (4) 

G {a} \<ij> i J 

where the connectivity matrix of the lattice is Gy. The dynamics of the lattice, represented by the 
sum over random (|)^^ graphs or their dual triangulations X^c takes place on the same timescale as the 
dynamics of the spins, which is manifest in the numerical simulations where the local "flip" moves that 
change the lattice connectivity to perform the sum ^^ are mixed in with spin cluster updates. An 
interesting related question, to which the answer is not immediately clear, is what exponents would be 
measured in a simulation on a fixed 2D quantum gravity graph. 

A first refiex is to guess that with a fixed lattice we would find the Onsager exponents. This appears 
to be correct when we consider random lattices with "weak" disorder - ie sufficiently close to regular 2D 
lattices, where the Harris criterion (roughly paraphrased as "weak disorder has no effect if a < 0") gives 
the Ising model as a marginal case, and simulations come down on the side of Onsager exponents [p[ n. 
Other fixed random lattices which are close in some sense to flat lattices Q, such as Voronoi lattices also 

^The topology of the underlying lattice is not expected to affect the spin model critical exponents, though it does have 
an effect on the exponent "y string- _ 

^The analytical work of the Dotsenkos |9|, however, suggests that randomness in the bond values can weaken the specific 
heat singularity. 

^Having a fractal dimension of 2, for instance. 



give Onsager exponents ||l^. There has been one simulation to date on a single random (fP' graph [Q 
generated by the Tutte algorithm which proved an efficient means of producing large pure (ie no matter) 
2D gravity lattices [|2| , that foimd results that appeared to be compatible with the Onsager exponents. 
A moment's reflection suggests that this is a rather surprising result - the fractal dimension of 2D gravity 
coupled to conformal matter with central charge c has been calculated |I3], and modulo some caveats 
about lattice dependent effects on the lower moments of the density of nodes ll4] , is given by 



-H>.^)f (^) 



where Q = ^(25 - c)/3 and 7 = (Q - ^Q"^ - 8)/2. We find that D for c = 1/2, the fsing model 
interacting with a dynamical lattice, is equal to 3 and D for c = 0, ie pure 2D gravity, is (5 + vl3)/3 — 
2.8685... We would thus expect to find d ~ 2.8685 for the graph used in [0|, which is very different from 
the 2 of a flat lattice. Unfortunately, we can't deal analytically with a single fixed (fr' graph to see what 
exponents we might expect, but we can aspire to taking a quenched average over a set of different fixed 
(j)^ graphs. 

It is often the case in statistical mechanics that quenched averages, in which the disorder is frozen 
in on the timescale of the spin dynamics, are more interesting than annealed averages, particularly in 
systems where frustration is present such as spin-glasses fl^. All of the couplings are positive in the 
model of equ. M) , so we have no frustration and hence would not expect to see any of the rich behaviour 
of [ p5[ , but it is still interesting to ask what happens to the critical exponents of the model when we 
consider a quenched sum over the random graphs rather than an annealed sum. In a simulation we 
would implement such a quenched summation by simulating our Ising model on an appropriate set of 
fixed random graphs and then taking a cross lattice average to obtain our observables, as is done in the 
simulation of spin glasses to average over different sets of spin couplings. In ||l^ it was suggested that 
applying the replica trick |15|] in conjunction with the continuum methods of [y allowed the calculation of 
the exponents in a quenched average over 2D gravity graphs. For a quenched summation the free energy 
F is calculated on each graph before summing over the random graphs 

F = Y,\ogZ{G), (6) 

G 

where Z(G) is the partition function for the Ising model on a particular random graph 

Z{G) = E e^P /5 ( E ^^^■^''^^ + H E ^^ ] • (7) 

{ct} \<u> » / 

The awkward logarithm inside the sum can be dealt with by replacing \ogZ ~ [Z"" — l)/n, n ^ to give 
a more familiar partition function, now with n replicas of the original matter. Observables are calculated 
with this replica partition function and the limit n ^ is taken in order to obtain the desired quenched 
averages. Translating this approach directly to the continuum language of [y , we can obtain the quenched 
value of the one point function i^$(A) of an operator <&, the magnetization or the energy for instance, 
by inserting one copy $° of $ into the fixed intrinsic area partition function Zn{A) for n copies of the 
matter theory 



/ D(j)DX exp(-S'")5 (/ e^'f'y/gd'^z - A) J '^'^ eft' ^/gd'^ z 



F^{A) 



Zn{A) 

n^Q (8) 

where the action S*" is composed of a Liouville part and a Feigen-Fuchs action |1^ for n matter copies 

^M^^iZJ "^'^ (dX'^dX'^ + '-ao^RX^^ (9) 



a=l' 



and uq is determined by c = 1 — 12q!q. Applying the standard scaling arguments of DDK j|] to such 
a theory and taking the limit n — > to get the quenched averages, we find the following values for the 
quenched conformal weights 



^quenched 



±Vl + 24Ao-l 



(10) 



If we plug in the values for the energy and spin operators in the Ising model we find the following 
quenched weights 



A„ 
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0.6513878 



0.1452847, 



fll) 



which, using equ.(0), gives the following values for the critical exponents a and (3 

6-2^/l3 



/? 




-0.8685170 
0.4167516. 



fl2) 



Applying the scaling relations then gives the full set of quenched critical exponents, which we repeat from 
|lq| in Table. 1 below for completeness, along with the fixed lattice (Onsager) and annealed (dynamical 
lattice KPZ/DDK) exponents. 
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Quenched 


-0.8685169 


0.4167516 


2.0350137 


5.8830375 
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2 - 0.7094306d 



Table 1: Critical exponents for Ising models 

We thus predict a set of quenched exponents that are close, but not equal to, the dynamical lattice 
exponents. The non-rational exponents might be a little disquieting at first site, but they seem to be a 
consequence of forcing the model to live in a gravitational background with the "wrong" value of c (ie 
c = effectively as we have a total central charge of nc, n — > in the quenched theory rather than the 
natural c = 1/2). Similar non-rational exponents were calculated and observed for percolation, where the 
natural value of c is eaual to zero, on c 7^ string worldsheets in p8[ . In fact if we insert the appropriate 
values for the fractal dimension calculated from equ.(P) into the above table we see that we obtain the 
same values (1, respectively) for the exponents v and 77 in both the annealed and quenched cases. If 
we take these two exponents as our starting point and calculate the rest from the scaling relations, we 
could argue that the quenched model differs only from the annealed model in living in a slightly different 
dimension. We have not found the fixed lattice values for the exponents in the above calculation because 
switching off 2D gravity completely corresponds to the limit c -^ —00 rather than c — > 0. The effects of 
2D gravity are thus still visible, even in an ensemble of fixed <^^ graphs. 

So much for the theory. The results of |ll|] suggest that the Onsager exponents appear on a fixed <^^ 
graph, so an average over a set of such exponents would still give Onsager exponents, in spite of the above 
arguments. In this paper we report on a simulation where we implement a quenched average explicitly by 
simulating O(IOO) Ising models on different fixed <^^ graphs simultaneously using a parallel computer. We 
have resurrected and modified the program used in earlier simulations on dynamical lattices [^ for this 
purpose, so we perform a microcanonical (fixed number of nodes) simulation on 0^ graphs of spherical 
topology without tadpoles or self-energy bubbles with iV — 124, 250, 500, 1000 and 2000 nodes. We used 
the Tutte algorithm to generate the one hundred or so (\? graphs of each size that were required. We 



simulated at a range of /? values from 0.1 to 1.1, using hot and cold starts where appropriate, with lOK 
metropolis thermalization sweeps and 50K Wolff cluster algorithm update sweeps at each (3 value and 
graph size. Note that the absence of NFLIP = N dynamical flip moves between each spin update in 
the quenched simulation means that our statistics are rather more modest than the dynamical lattice 
statistics in p| . The various observables were measured on each graph and a cross graph average taken to 
implement the quenched sum. The errors quoted are those coming from the average across the different 
graphs, rather than individual graph statistical errors. 

The standard observables were measured: the energy E, specific heat C, magnetization M, suscepti- 
bility X) correlation length ^ and Binder's cumulant for the magnetization 

along with the acceptance for the Wolff algorithm, to verify that an appropriate number of sweeps was 
being carried out. The number of Wolff sweeps was adjusted so that (1 / average size of Wolff cluster) 
sweeps were performed per measurement at each (3. In our original dynamical lattice simulations in 
P] we had some difficulty in finding a reliable estimate for the crossing point of Binder's cumulant for 
different lattice sizes, which is the usual means of locating the critical point with these methods, so we 
adopted an alternative approach that gave good results for the exactly soluble case of the Ising model, 
where the critical point is known analytically. The expected scaling behaviour for Binder's cumulant is 
U ~ U{tN^^'^^), so its derivative with respect to t, or equivalently /?, will scale as l/vD at [3c- We would 
also expect the maximum slope to scale in the same way 

maxf— ^ ~ TVi/''^ (14) 

\dl3 J 

which can be used to obtain uD. Knowing this the finite size scaling relation 

\I3^ -I3^\ ~ TV-i/"^^ (15) 

with (3-^ obtained from either the peak in the specific heat or the maximum slope in the cumulant, can 
be used to extract /3^. In Fig. 1 we show the maximum slope of the cumulant plotted against TV on a 
logarithmic scale along with the fitted line with slope \/vD — 0.32(1). We thus have vD ~ 3.1(1), which 



can then be fed into the finite size scaling relation in equ.(15) to get (3^ . In Fig. 2 we show the results 
for both the maximum slope in the cumulant and the peak in the specific heat, giving /3^ = 0.778(2), 
where we have dropped the smallest lattice size from the fit to the cumulant slope. We also found that 
it was possible to extract 13c from the crossing of cumulants for different lattice sizes with the current 
data, which gives an additional estimate of (3c ~ 0.78(1) compatible with the value above, as is the 
corresponding estimate of vD with this method. 

The scaling relation a = 2 — vD thus gives a = —1.1(1) for the specific heat critical exponent, which 
is a long way from the logarithmic scaling [a = 0) given by the Onsager exponents. It is also possible to 
attempt to extract a directly from the singular behaviour of the specific heat 

C ~ B + Cot-" (16) 

or from finite size scaling 

C ~ B' + C'o N"/"". (17) 

Using the singular behaviour proves to be hopeless with our data as essentially any value for a can be 
obtained with comparable x^s, but the finite size scaling fit gives a minimum x^"fit for B' = 1.8 and 
a/vD = —0.37(1), which gives a — —1.1(1), in good agreement with the value obtained indirectly from 
the scaling relation. 

A measurement of the exponent 7 governing the divergence of the susceptibility is less useful in 
distinguishing between the different sets of exponents as the quenched and annealed values are essentially 
identical (~ 2) and the Onsager value is close to these (1.75). A fit to 

X - Xo t--^ (18) 



gives 7 = 2.1(2) on the N = 2000 lattices and 7 = 2.0(2) for N = 1000, which although it favours the 
dynamical/quenched values, is not really sufhcient to exclude the Onsager values. The finite size scaling 
fit gives a much poorer fit with an even larger value of 7. In Fig. 3 we plot the divergence in x for the 
two lattice sizes along with lines of slope 1.75 and 2. It is clear from the diagram that it is indeed not 
really possible to decide conclusively between the candidate values from the data. Fitting (3 is also rather 
unilluminating, whether by using M ~ Mq t^ or M ~ Mg N~^/'^^ as we find values of 0.2-0.3 stranded 
in no man's land between Onsager and dynamical/quenched exponents. We found in [Bl that the fits to 
/? on dynamical lattices were the poorest, so a similar result here is not surprising. It is not clear why 
the fit to the magnetic exponent should be so poor, but it was pointed out in H] that the dynamical 
critical exponent z for the magnetization was much larger than that for the energy on dynamical lattices 
(with the Wolff algorithm) possibly because of clusters becoming trapped in "bottlenecks" in the graph 
geometry. It is therefore conceivable that one has to work harder to overcome autocorrelation and finite 
size effects in the measurements of /? than for other exponents, whether on dynamical or quenched 2D 
gravity graphs. 

We can also fit the exponent i' by looking at the divergence of the correlation length ^ at /3c. As this 
involves a double fit, first for the correlation length 

C\j = < ^^X^^»crj'5(f^y - r) > ~ cxp(-^) (19) 

(where dij is the geodesic distance on the lattice between points i and j and n{r) = ^^ S{dij — r)) 
and then to the divergence in ^, the fit is not likely to be very accurate. Doing this however, we find 
v = 0.87(2) for the data on the N = 2000 graphs, which is identical to value of v we measured on 
dynamical graphs in g. If we take the theoretical values for D for both the quenched and annealed 
graphs we would expect v = 1, just as for the Onsager exponents. 

In summary, the fits described above are thus compatible with both those in ||5[] on dynamical (p'^ 
graphs and with the similar quenched exponents. The best results are for i^D and a/vD, both of which 
give a cusp (a ~ —1.1(1)) in the specific heat, rather than a logarithmic divergence, as expected for a 
flat 2D lattice. The measured value of 7 on the largest lattices used, 2.1, is close to the ~ 2 expected 
on dynamical/quenched graphs, but the Onsager value of 1.75 lies not far outside the error bars. The 
measured value of v is identical to that found for dynamical graphs, but one cannot use v to distinguish 
between Onsager and dynamical exponents. The values of (3 measured are as disastrous as those found in 
|q] on dynamical lattices, and hence not very useful. Finally, it is worth noting that the measured value 
of /3^ in the quenched simulation is close to the calculable (3c on dynamical graphs (0.7733185). 

Although it would be dangerous to claim that there was conclusive evidence for non-Onsager exponents 
in the quenched ensemble on the basis of fits to vD and a alone, the data for these is at least consistent 
with this, and none of the other fits contradict it. Unfortunately, the results are nowhere near accurate 
enough to distinguish between the postulated quenched exponents and dynamical graph exponents, even 
with the 0(100) different graphs used. If we accept that the quenched ensemble does give non-Onsager 
exponents we are left with the task of explaining the results of |0] on a single 2D gravity graph, where 
acceptable fits were found to the Onsager values. It is possible that taking the continuum limit in two 
different ways, on a fixed graph and in an ensemble of fixed graphs, simply gives different exponents, 
though as one expects self-averaging for quantities such as the energy and magnetization it is difficult 
to see how this might happen. Another possibility is that the particular random graph chosen in pTl 
produced a set of pseudo-Onsager exponents by chance, as there appears to be a wide variation in the 
apparent critical points and exponents if these are fitted in the individual graphs in our ensemble. Earlier 
simulations on fixed random lattices, of the XY model in particular |19| , have shown that finite size 
effects can be particularly important too. In closing it is amusing to note that, formally at any rate, 
we would expect to see the same quenched exponents appear in a simulation of an Ising model on a 
single dynamical graph in which the back reaction of the spins on the lattice was switched off (by always 
carrying out flips rather than performing a Metropolis test on the resulting energy change) as this also 
has c = 0. 
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Figure Captions 

Fig. 1. The maximum in the slope of Binder's cumulant is plotted against the lattice size to extract lyD. 
The best fit is shown as a line. 

Fig. 2. The (3 value of the maximum in the slope of Binder's cumulant and the peak in the specific heat 
are plotted against N^^^'^^ to extract the critical coupling /3^, given by the intersection of the 
best fit lines with the y axis. 

Fig. 3. A finite size scaling plot of x, with lines of slope 1.75 (the steeper) and 2 added for comparison. 
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